Numerical analysis of the quantum dots on off-normal incidence ion sputtered surfaces 
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We implement substrate rotation in a 2+1 dimensional solid-on-solid model of ion beam sputtering 
of solid surfaces. With this extension of the model, we study the effect of concurrent rotation, as the 
surface is sputtered, on possible topographic regions of surface patterns. In particular we perform a 
detailed numerical analysis of the time evolution of dots obtained from our Monte Carlo simulations 
at off-normal-incidence sputter erosion. We found the same power-law scaling exponents of the 
dot characteristics for two different sets of ion-material combinations, without and with substrate 
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I. INTRODUCTION 

The size-tunable atomic-like properties of (e.g. II- VI 
and III-V) semiconductor nanocrystals have diverse ap- 
plications. These properties arise from the quantum con- 
finement of electrons or holes in the quantum dots (QDs) 
to a region on the order of the electrons' de Broglie wave- 
length. Examples of the applications can be found in 
solid-state quantum computation^ and in electronic and 
opto-electronic devices like diode lasers, amplifiers, bi- 
ological sensors, electrolumniscent displays, and photo- 
voltaic cells .-^ 

A more recent method of fabrication is the sputtering 
of semiconductor surfaces with a beam of energetic ions 
impinging at an angle 9 with respect to the direction 
perpendicular to the surface.'^ This has been shown to 
be a cost-effective and more efficient means of producing 
uniform high-density semiconductor nanocrystals,— i^i^ 
in contrast to previous methods such as epitaxy, litho- 
graphic techniques, colloidal synthesis, electrochemical 
techniques, and pyrolytic synthesis. 

Using the continuum theory, it was shown that QD 
formation by = sputtering is restricted to a very nar- 
row region of the parameter space.-'? It has also been 
shown that dot formation is possible for > 0, within a 
broader region, under concurrent sputtering and sample 
rotation.^ Using a simple discrete solid-on-solid model, 
which includes the competing processes of surface rough- 
ening via sputtering and surface relaxation via thermal 
diffusion, we recently found that for > 0, without sam- 
ple rotation, a dot topography is only one among other 
possible topographies which may arise. The type of the 
emerging topography depends on the longitudinal and 
lateral straggle of the impinging ion as it dissipates its 
kinetic energy via collision cascades with atoms within 
the material.— 

In this study, we implement sample rotation in the sim- 
ulation model along the lines of Refs. 11. and 12. For vary- 
ing values of the ion parameters we find different kinds 
of surface patters, including dots. The dots are similar 
to those obtained without sample rotation, but the un- 
derlying oriented ripple structures are lacking. We study 
the time evolution of dots emerging from oblique ion in- 



cidence in more detail, performing simulations for two 
different sets of parameters, corresponding to two differ- 
ent materials (GaAs sputtered with Ar and Si sputtered 
with Ne) with and without rotation. We found that with- 
out rotation the average number of dots decreases with 
increasing fluence, but stays approximately constant for 
a rotated sample. Furthermore the uniformity of dots is 
greatly enhanced by sample rotation. Remarkably, both 
materials exhibit the same scaling of the dot character- 
istics with sputter time. 

The rest of the paper is organized as follows: In section 
im we shall briefly review the continuum theory of sput- 
tered amorphous surfaces. In section IIIII we shall give 
a brief description of the discrete simulation models ap- 
plied in this work. In section IIVI we study the effect of 
rotation on the possible topographies reported in Ref. [lol : 
expecially off-normal incidence dot formation. Finally, in 
section IVl we shall present and discuss our results of the 
analysis of dot topographies. 



II. CONTINUUM THEORY 

In the seminal theory of P. Sigmund on ion-beam 
sputtering of amorphous and poly-crystalline targets, 
it was shown that the spatial energy distribution -E(r) 
of an impinging ion may be approximated by a two- 
dimensional Gaussian of widths a and /i, parallel (z'-axis) 
and perpendicular (a;', y'-axes) to the ion beam direction, 
respectivel y -^^d'* , 



■ exp 



X'2 + y'2 



, (1) 



where e is the total energy of an impinging ion and a the 
average penetration depth. 

When describing the surface by a two-dimensional 
height field h{r,t) (solid-on-solid description), the nor- 
mal erosion velocity v = -[l + {'Vh)'^]-^^'^dth at the point 
r = (0,0,0) is proportional to the total power trans- 
ported to this point by the ions of the impinging ion 
beam. It may be expressed a o^^d^ 



/ dr$(r)^(r). 
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The integral is taken over the region R containing all 
the points at which the deposited energy contributes 
to the total power at r = and pc is a proportional- 
ity constant r^*' <I>(r) is a local correction to the uniform 
flux Jf. Note that shadowing effects among neighboring 
points and redeposition of eroded material are ignored. 

By standard arguments, an equation for the evolution 
of the surface height field h{x, y, t) due to sputtering can 
be derived from Eq. ([2]), which takes on the forniiiii^ 

{dth)s ~vf + v'fdxh + Vxdlh + Vyd^h 

+ ^{d,hf + ^{dyhf +ils{v.t), (3) 

where the x-axis is parallel to the ion beam direction, 
?;/ is the erosion velocity of a flat surface and Vx{9, u, /x) 
(T, /i)] is the surface-tension coefficient along [per- 
pendicular to] the ion beam direction. Depending on the 
parameters 9, a, fi, the quantity can exhibit positive 
as well as negative values, whereas Vy is always negative. 
Xx and Xy are the coupling constants of the dominant 
non-linearities along the respective directions, rjs is an 
uncorrelated noise term, with zero mean. This term rep- 
resents the random arrival of the ions onto the surface. 

The surface height also evolves due to surface particles 
hopping from one point to the other. On a coarse-grained 
level this can be described by the continuity equation for 
the conserved particle current 

{dth)D = -y ■i+VD{r,t), (4) 

where the local current density j is a function of the 
derivatives of h; i.e, j = j[V'"ft., (Vft.)"] (m and n are in- 
tegers), and the acceptable functional form is subject to 
symmetry constraints. rjD reflects the randomness inher- 
ent in the surface-diffusion process. A commonly stud- 
ied example of a surface migration model is the MuUins- 
Herring model, which leads to the particularly simple 
form V • j = Z?V^/i for the current density. 

Thus, considering Eqs. ^ and (|3]), the time evolu- 
tion of a sputtered surface is governed by a Kuramoto- 
Sivashinsky type equation 

dth = ~vf + v'fdxh + v^dlh + vydlh + ^{dxhf 

+ ^{dyhf ~DV^h + r,{v,t). (5) 

The exact form of each coefficient of Eq. [5] is given in 
Ref. [l^ and will be used later in the analysis of our results 
(see Fig. (P in section dllH '). 

The anisotropy {v^ ^ Vy, Xx 7^ Xy) arises from the 
oblique incidence which leads to different erosion rates 
(or different rates of maximizing the exposed area) along 
parallel and perpendicular directions relative to the ion 
beam direction. From the linearized Eq. ([5]) it is easy 
to see that periodic ripples are growing with wavevector 
k ■^max(|i^^/y)|/2Z) driven by the dominant negative 
surface tension. These ripples are stabilized by the sur- 
face diffusion term. With prolonged sputtering surface 



slopes become too big to neglect the nonlinear terms and 
the ripples are destroyed, after which a new rotated rip- 
ple structure may emergeJ^ 

If the sputtered substrate is rotated with sufficiently 
large angular velocity, or if = 0, no externally induced 
relevant anisotropy is present. In this case ripples do 
not develop for the class of materials we study here and 
the relevant evolution equation is the isotropic version of 
Eq. ([5]) , which predicts cellular structures of mean width 
~ {Dl\v\y'^'^ and mean height ~ I'^l/lxl'^ According to 
Ref.[^ these cellular structures eventually evolve to either 
a dot topography (x >0) or a hole topography (x <0); 
which explains the normal incidence dot^^S^Z, and hol o^°i^^ 
topographies found in previous experiments. 

Thus, one would expect oblique incidence dot forma- 
tion to result from isotropic sputtering induced by sub- 
strate rotation. However, in our recent simulationsiS we 
found dot topographies without sample rotation, though 
with some underlying (anisotropic) ripples oriented par- 
allel to the ion beam direction. We found that the oblique 
incidence dots arise from a dominance of the ripples ori- 
ented parallel to the ion beam direction, over those with 
a perpendicular orientation; the dots being the remains 
of such perpendicular ripples (see Fig.[2]below) after long 
times. 

In the next chapter we will outline the simulation 
method we have used in Ref. [13 and which we have ex- 
tended to the case with rotation as presented in this work. 



III. MONTE CARLO SIMULATION 

To simulate the competing processes discussed in sec- 
tion [n] on a discrete system, we use the Monte Carlo 
model of sputter erosion introduced in Ref. [13 (HKGK 
model). To make the paper self-contained, we provide 
some details here. We simulate the sputtering process on 
an initially ffat surface of size with periodic boundary 
conditions, by starting an ion at a random position in 
a plane parallel to the initial surface, and projecting it 
along a straight trajectory inclined at angle 9 to the direc- 
tion perpendicular to the averaged surface configuration 
and at an azimuthal angle (p. An ion penetrates the solid 
through a depth a and releases its energy according to 
Eq. IH). An atom at a position {x,y,h) is eroded (see 
Fig. [T|) with probability proportional to E{r') given in 
Eq. 11]). It should be noted that, in accordance with the 
assumptions of the theoretical model o^'^i^^'^^ , this sput- 
tering model neglects evaporation, redeposition of eroded 
material and preferential sputtering of surface material 
at the point of penetration. The surface is defined by a 
single valued, discrete time dependent SOS height func- 
tion h{r, t) — h(x, y, t). The time t is measured in terms 
of the ion fluence; i.e, the number of incident ions per 
two-dimensional lattice site {x,y). The choice of the pa- 
rameters a, fi and 9 is discussed below. Substrate ro- 
tation, at constant angular velocity u > /D, is imple- 
mented by keeping the substrate static and then rotating 
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the ion beam instead, which is equivalent to keeping the 
ion beam direction fixed and rotating the substrateii. In 
our simulations we choose (j) randomly from the interval 
< < 27r, since, in the large uj limit the solid appears 
to be sputtered from all angles (j)^ 
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FIG. 1; The model consists of a square field of discrete height 
variables h{x,y), corresponding to piles of h{x,y) particles 
at position {x,y). Here a projection to the h — x is shown. 
Left: Each ion impact is modeled by an distribution describing 
the energy deposited by the ion. Atoms on the surface are 
removed with a probability proportional to the energy. Right: 
Surface diffusion; by decreasing height differences the energy 
is decreased. 

Our model of the sputtering mechanism sets the time 
scale of the simulation in a way, which allows direct com- 
parison with experiments. Any relevant surface diffu- 
sion mechanism may be combined with this sputtering 
model. Here, we use a realistic solid-on-solid model of 
thermally activated surface diffusion^"^. Surface diffusion 
is simulated as a nearest neighbor hopping process with 
an Arrhenius hopping rate 



R{E, T) = RQe:cp{~E/kBT), 



(6) 



where i?o = ksT /h^ and T is the effective substrate tem- 
perature (see below). The energy barrier E = Eyn + 
n-nEin + Ese consists of a substrate term (Eyn =0.75 
eV), a nearest neighbor term (£';„ =0.18 eV), and a step 
barrier term {Ege =0.15 eV), which only contributes in 
the vicinity of a step edge, and is zero otherwise. In 
each diffusion sweep all surface particles are considered, 
and those that are not fully coordinated may hop to a 
neighboring site. Note that we have to use a higher ef- 
fective temperature^! in our simulation in order to ac- 
count for the greatly enhanced surface diffusion due to 
thermal spikes. A thermal spike is a series of sharp peaks 
and minima in the spatio-temporal distribution of the 
surface temperature, arising from the occurence of lo- 
cal heating of the surface right after every ion impact, 
followed by rapid cooling. Hence, we have used a higher 
effective temperature fc^T =0.1 eV, which was estimated 
in our previous work and which roughly corresponds to 
experimental sputtering at room temperature. 

The model captures the essential features of sputtered 
surfaces at nanometer lengthscales; especially nonlinear 
effects^i^i. Although a direct mapping of this model 
to the continuum equations is unavailable, the model 



is expected to be consistent with variants of the KS 
equation ] 



IV. EFFECT OF ROTATION ON THE 
TOPOGRAPHIES 

In this study, a lattice of linear size L = 128 is used 
with a lattice spacing corresponding roughly to a dis- 
tance of 0.5 nm. The simulated time is set such that 
1 ion/atom corresponds to an ion fluence of 3x10^** 
ions/cm^. We use a sputter yield of about 7 surface 
atoms/ion, as compared to 5 Si02 molecules/ion (re- 
duced to 0.1 Si02 molecules/ion with H ion) and 0.3-0.5 
molecules/ion in the experiments of Refs. [2y and [27l re- 
spectively; this may result in lengthscales differing from 
those of the cited experiments, but we found in our previ- 
ous studieaiSi22ii^ that predicted universal features exist, 
which are directly comparable to experimental results. 
We will subsequently discuss such features for rotated 
samples. 

In Ref. [13 six regions with different topographies (in- 
cluding e.g, smooth, hole, and dot regions) were shown 
to emerge for ion collisional parameters a = 6, < cr < 5, 
< /i < 5 (all measured in units of lattice spacings), at 
time t = 3 ions / atom. In this section we study the effects 
of sample rotation on these topographies, in particular ef- 
fects on the region with dot topography. The differences 
in topographies between rotated and unrotated samples 
are visualized in Fig. [2] and Fig. [3l for a = 6, i = 3, and 
9 = 50°. 

As seen from Fig. |3] (and Fig. [H see below), no 
anisotropy can be found with substrate rotation, as ex- 
pected from the continuum theory. The ripple structures 
obtained for /i < 2 (Fig. [5]) do not appear for rotated sub- 
strates. The underlying parallel ripples of the dot region 
(topmost row of Fig. [2]) are also absent for rotated sub- 
strates. However, hole formation is not suppressed, we 
get holes with as without rotation as is visible in Fig. [21 
This fact can be understood from the continuum theory, 
which predicts roughly equal erosion rates along both di- 
rections for parameters in the hole region, hence there 
is no anisotropy to be destroyed. Furthermore, ripple 
patterns perpendicular with resepct to the ion beam di- 
rection are replaced by non-oriented structures, and the 
ordered parallel ripples are no longer present if the sub- 
strate is rotated, (see Fig. 13]), 

For a closer inspection, we calculate the structure 
factors, 5'(k, t) = |/i(k, f)p from the Fourier transform 
/i(k, i) of the height field h{r,t). In particular we con- 
sider four prototypical topographies marked by letters S, 
H, N, D in Fig. [31 S stands for "relatively smooth" , H 
for "hole", N for "non-oriented structures", and D for 
"dots" . The results are shown in Fig. [H As can be seen 
from this figure, and as expected, there is no anisotropy 
visible in all cases. In the case of the relatively smooth 
surface S, there is also no characteristic lengthscale. For 
the hole topography, H, there is still no specific length- 




13 5 



a — I- 

FIG. 2: Possible topographies of the model, within the ex- 
perimental constraints considered in Ref. [lO- t — 3, a — Q, 
6 = 50°. Left - right columns: a = 1, 3, and 5, respectively. 
Bottom row - top row: fi — 0.5, 1, 1.5, 2, and 5, respectively. 
The last two profiles of the top row belong to the dot region 
(region V) of Ref. M- 



scale but there now exists an upper bound kub on |k| due 
to the presence of the holes. On the surface with non- 
oriented structures (N) a well defined lengthscale with 
kub as well as a lower bound kib can be found. And fi- 
nally, in the case of the dot topography (D), we also have 
a characteristic lengthscale, but kib is shorter here than 
for the N topography, which implies that the average sep- 
aration of the dots is larger than that of the non-oriented 
structures, as expected from Fig. [3] 

According to the continuum theory, there exists a sin- 
gle effective surface tension coefficient Vav = Vx + Vy and 
a single nonlinear couplingi^ Xav for all directions in ro- 
tated samples, since there is no anisotropy left in the sys- 
tem. Fig. [5] shows a plot of these coefficients based upon 
the explicit expressions in Ref. [l6|. As can be seen from 
Fig. [5] and Eq. ([5]) the surface roughens with time, with 
smaller Vav {i^av <0) corresponding to higher roughness. 
Surfaces in the parameter range for which 0> Vav S> -1, 
are relatively smooth (S) topographies. 
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FIG. 3: Profiles obtained from simultaneous sputtering and 
rotation, using the same parameters as in Fig. [2] t = 3, a = 
Q, 9 = 50°. Left - right columns: a = 1, 3, and 5, respectively. 
Bottom row - top row: ^ = 0.5, 1, 1.5, 2, and 5, respectively. 
The last two profiles of the top row belong to the dot region 
(region V) of Ref. [l^. Structure factors of the lettered profiles 
are provided in Fig. H] [S - (relatively) smooth; H - hole; N - 
non-oriented structures; D - dot]. 




S H N D 

FIG. 4: Structure factor of the lettered surface profiles in 
Fig-E] S (relatively) smooth; H hole; N => non-oriented 
structures; D dot. k = at the centre; and k — ^(±1, ±1) 
at the corners. 



Considering Fig. [5l one sees that \vav\ first increases, 
in accordance with the increasing height difference on 
the greyscale charts on the profiles of Fig. [3] and then 
it decreases as we tend towards the dot region. These 



FIG. 5: The coefficients, Vav and x^v, of the isotropic version 
of Eq. [5l^^ for the rotated case, as functions of /i. (a) a = 
1.0, (b) a = 3.0, and (c) a = 5.0. 




FIG. 6: Surface roughness, W, with dots excluded. Main 
plot: W as a function of /i for a — 5. Inset: W as a function 
of a for n — 5. t — 3 ions/atom. 



changes in the roughness are not visible on the greyscale 
charts due to the appearance of the dots, which are con- 
siderably higher than an average surface protrusion in 
the dot-free profiles. Therefore we have also studied the 
surface roughness, with the dots excluded, as shown in 
Fig. [6] (Details of our dot-isolation method are discussed 
in the next section). In this figure, the roughness W as 
a function of (cr — 5) is shown in the main plot, where 
the roughness first increases and than decreases again, in 
accordance with the continuum theory. The inset shows 
a plot of W versus cr, for /i = 5. 

When the local surface slopes become significant (with 
prolonged sputtering), nonlinearities become relevant. It 
has been shown that for 6* = crossover to the nonlinear 
regime either gives rise to dot formation (if Xav > 0) or 
hole formation (in case Xav < 0) if ion-induced effective 
surface diffusion is the dominant relaxation mechanism^. 
This is consistent with our results for parameters at 
which Xav > in Fig. [5] (i.e, for /i > 3). We also 
found holes for Xav < ("H" region), but hole forma- 
tion for long times is not as widespread as Fig. [5] seems 
to indicate. In particular, the hole topography eventually 
evolves into cellular structures similar to those shown in 
Fig. [7] at long times. Since Vav 7^ 0, the surface rough- 
ening is not wavelength independent, which explains the 
presence of the non-oriented protrusions. 

For very small longitudinal and lateral straggle, cr < 
1, /i < 0.5, i.e. in the "S" region, we did not find any 




FIG. 7: Time evolution of the relatively smooth topography 
of Fig. [2] with rotation, a = 3, fj, = 0.5. (a)-(d): t = 3, 40, 
90, and 150, respectively. 



structure up to the longest simulation times. Note that 
with increasing a, the l^avl ~ interval is reduced. 

For non-oriented structures "N" , simulations at longer 
times reveal only slight changes in the structures; no dot 
formation (see Fig. [5]) appears. Hence, we have only ob- 
served dots with rotation wherever they can be found 
without rotation. 




FIG. 8: Time evolution of the non-oriented structures arising 
from rotation of rippled region of Fig. (2] a = 3, /i = 1.5. (a) 
-(d): t — S, 40, 90, and 150, respectively. 



When considering the emerging topographies of Fig. [3] 
at other angles of incidence, we found notable changes 
with 6, as illustrated in Fig. [9] for cr = 1, /i = 2 and 
a ~ 1, fi = 5). That is, changes from smooth — > hole 
non-oriented structure, and back to smooth (no struc- 
ture) topographies appear. These changes imply that the 
roughness initially increases and then starts to decrease 
with increasing 9 as shown in Fig. 1101 In this figure 
cr = 1 and data for /x = 2 and /it = 5 are represented 
by (black) circle and (red) square symbols, respectively. 
The behavior of the roughness with varying 9 shown in 
Fig. [TU] is in agreement with the experiments reported 
for Ar-ion sputtered rotating InP surfaces in Ref. 5j. But 
note that the roughness data reported in this experiment^ 
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"in lattice units 



FIG. 9: Topography changes as a result of changing the angle 
of incidence 6, for a = 1; ^ = 2 (top row), and 5 (bottom 
row). Top row, L-R: 9 = 10, 30, 40, 50, 80. Bottom row, L-R: 
e = 5, 10, 30, 50, 80. 



were obtained for an exposure time of 1200 sees., i.e. in 
the steady state regime. To use such steady state val- 
ues in Fig. [TO] would require simulation times beyond our 
computational time constraints. 




FIG. 10: (Color online) Surface roughness as a function 
of 9 for profiles shown in Fig. [9] Data for fi = 2 and 5 are 
represented by circle and square symbols, respectively. In 
both cases a — 1, t = 3 ions/atom. Error bars are included, 
but they are much smaller than symbol size except at = 10° 
{fj, = 5). For this case for about half of the runs the structure 
had already evolved to a non-oriented structure pattern with 
a considerable roughness, while the other half exhibited still 
a rather smooth surface. Hence, there seems to be a sharp 
transition in time from smooth to rough, where the transition 
time of finite samples fluctuates strongly. 



V. NUMERICAL ANALYSIS OF THE DOTS 

Focussing on topographic region "D" where dots have 
been found also without rotation, we obtained the ex- 
perimental parameters corresponding to this region from 
SRIM simulations^^ as shown in Table HI We have per- 
formed two sets of simulations using the SRIM results: 
(i) 1.6 keV Ar ion sputtering of GaAs surfaces; where a 
= 6, a = 5, fi = 3.8. (ii) 650 eV Ne ion sputtering of 
Silicon surfaces; where a = 6, a = 4.4, /i = 3.2. In both 



cases we have used an ion beam inclination of 6 —50° to 
the vertical. 

In order to study the time evolution of dot character- 
istics (e.g. dot density, area, and height), we have used 
a similar clustering approach as in Ref. |24| . A dot is de- 
fined as a cluster of points of local height maxima. This 
is done in two steps 

• For any given time t (we do not mention the time- 
dependence explicitly here) we consider the set of 
points 

M = {(x, y)\x, y e (1, • ■ • , L); h{x, y) > /i J, (7) 

where L is the linear size of the system, he is a 
cutoff height which the surface height at a point 
must equal/exceed for the point to be counted as 
(or part of) a doti^ We define he to be a function 
of the average surface protrusion, which has the 
form: he = hmin + p{{h) - hmin)- Where hrmn is 
the lowest surface height, (h) is the average surface 
height, and p is a fixed percentage. 

• We call two points in M neighbors if their distance 
is smaller than a given threshold dc, we use dc = 1 
here. Then, two points are called connected, if there 
exists a path from the first point to the other point 
such that all consecutive points along the path are 
neighbors. Now, each dot D is a subset of M (non- 
overlapping with any other dot) of maximum size 
such that all points in D are mutually connected. 
Hence, dots are the transitive closures of the neigh- 
bor relation on M. 

We start our simulations with a dot configuration ob- 
tained from a topography at < = 3 ions/atom and we 
choose a value p — po that yields the highest num- 
ber of sampled dots No (see Table |lT]). The initial 
clusters are shown in Fig. [Tl] for simulation of Ar-ion 
sputtering of GaAs without substrate rotation and with 
substrate rotation, respectively. The dot cross-section 
A{D) is defined A{D) = \D\; i.e the cardinality of D. 
The average height hd of a single dot D is defined as 
hdiD) = j^^J2{x,y)£D^x,y Notc that, in order to ex- 
clude non-dot surface protrusions (see top row of Fig. 
[TO)) from our analysis, we have included an upper bound- 
ary of 50 cluster points (i.e A{D) < 50) to our definition. 
In our analysis we only consider dots defined by these 
clusters. The results reported here are obtained from an 
average of 100 independent runs. 
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TABLE 11: Optimal parameters used to determine the cutoff 
height he 



sputtering of Si (Ne-Si). 



Simulation 


Po (%) 


No 


State 


Ar+ on GaAs 


50 


69 


unrotated 


Ar+ on GaAs 


10 


149 


rotated 


Ne+ on Si 


60 


76 


unrotated 


Ne+ on Si 


20 


139 


rotated 
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FIG. 11: (Color online) Sample surface profiles for the un- 
rotated case (top), and the rotated case (bottom) at t = 3 
ions/atom. In the figures on the right the clusters formed 
from the corresponding profile on the left (as defined in the 
text) are printed on top of the profile. The scales indicate the 
surface height measured from the lowest height. 



In Fig. [T^ we show sample surface profiles without 
rotation (top row), and with rotation (bottom row), for 
simulation of 1.6 keV Ar-ion sputtering of GaAs. As can 
be seen from the top row of this figure, the ripples (with 
parallel orientation to ion beam direction) that coexist 
with the dots become more ordered with time, whereas 
the dots decrease in number with time; more analysis of 
these dots is provided below. On the other hand (bot- 
tom row of Fig. [T^ . these ripples do not exist when the 
substrate is subjected to concurrent rotation, and the 
density of the dots is more uniform. 

In Fig. [T21 we present results of the average number of 
sampled dots Nc, and the average dot height He = hd, 
where the average is taken over all dots and all indepen- 
dent runs. The results of the average area of cross-section 
of the dots is presented in Fig. [Ml In Figs. [13] and [T4l 
open and closed circle symbols represent data obtained 
from Ar-GaAs with and without rotation respectively; 
while open and closed triangle symbols represent data 
obtained from Ne-Si with and without rotation respec- 
tively. The main result is that we found power law scal- 
ing of the dot characteristics with time; with or without 
rotation. In addition, we found the same scaling behav- 
ior for the two sets of simulation performed; i.e, 1.6 keV 
Ar ion sputtering of GaAs (Ar-GaAs), and 650 eV Ne ion 
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FIG. 12: Surface profiles for the non-rotated (top row) and 
the rotated (bottom row) case. L - R, t = 10, 20, and 30 
ions/atom respectively. The bar indicates the ion beam di- 
rection, and the scales indicate the surface height measured 
from the lowest height. 



To be more specific, as can be seen from Fig. [13] 
(a), the average number of sampled dots decreases with 
time as Nc ~ t^^ for non-rotated samples, where "0 — 
0.583±0.007, while it stays approximately constant for 
rotated samples ('0 ~ 0) . This indicates that the number 
of dots becomes insignificant without substrate rotation, 
so that one should use rotation if dot creation is the main 
purpose of the sputtering. Note that this result is already 
visible in Fig. [12] in qualitative form. 

On the other hand, as shown in Fig. [13] (b), the average 
height He increases with time as He ~ , the increase 
is more rapid for the unrotated case {g =0.409±0.004) 
than for the rotated case [g =0.159±0.007). As expected, 
the dot height is lower with sample rotation due to the 
enhanced smoothening effect of the rotation . This 
smoothening also explains the higher initial number of 
dots for the rotated case [Fig. [13] (a)]: The dots that are 
not "visible" in the unrotated case, since they do not 
surmount the ripples, become "visible" in the rotated 
case since the rotation prevents ripple formation. 

The mean cross sectional area Ac finally becomes al- 
most time independent for both cases. However, as 
shown in Fig. [14] , an initial power law scaling Ac ~ t~'' 
occurs, which is more significant for the rotated case 
(where = 0.113±0.006) than for the unrotated case 
(where <7 « 0). 

These results indicate that while no new dots are cre- 
ated with rotation (for > 0°), the uniformity of the 
density of the existing dots and the stability of the dot 
height are greatly enhanced with substrate rotation. Our 
results of the analysis of the dots are in agreement with 
previous experiments on dots obtained from normal inci- 
dence sputtering of semiconductors (Si),^ semiconductor 
alloys (GaSb, InSb);^'^ as well as oblique incidence dots 
obtained from simultaneous sputtering and sample ro- 
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observed experimentally so far. In particular we found 
20 



FIG. 13: (Color online) Time evolution of (a) the average 
number Nc, and, (b) the average height He (both in lattice 
units), of the sampled dots with time (in ions/atom). Sym- 
bols: open and closed circle (triangle) symbols denote data 
obtained with and without substrate rotation respectively, for 
Ar-GaAs (Ne-Si). 




FIG. 14: (Color online) Time evolution of the average area 
of crosssection Ac of the sampled dots with time for (a) the 
unrotated case, and (b) the rotated case. Symbols: open and 
closed circle (triangle) symbols denote data obtained with and 
without substrate rotation respectively, for Ar-GaAs (Ne-Si). 



tation (InP)j^ In these experiments, the dot height has 
been reported to increase with time; and the average dot 
size have been reported to become constant with time. 

To summarize, we have implemented substrate rota- 
tion for a solid-on-solid model of surface sputtering and 
used it to study the effects of concurrent rotation on the 
different possible topographies. In particular, we have 
studied the effect of rotation on the dot region as well 
as a detailed analysis of the time evolution of the dot 
characteristics (number, cross sectional area, and height) 
with and without substrate rotation. We found that dif- 
ferent materials whose sputtering parameters fall within 
this region exhibit the same scaling behavior. The num- 
ber of dots Nc formed in the absence of substrate rotation 
decrease with time as Nc ~ t~^-^^, whereas Nc is roughly 
constant with substrate rotation. Both with and without 
rotation the dot cross section A finally becomes indepen- 
dent of time, however, it initially decreases according to 
A with rotation. 

Additionally, for other choices of the sputtering con- 
ditions, we find different patterns which have not been 



transitions in time from one kind of surface structures 
(e.g. smooth, or holes) to other structures (like non- 
oriented structures), which can be explained only by the 
presence of non-linear effects. Hence, more sputtering ex- 
periments with different ion/substrate types and varying 
parameters are needed to verify whether the structures 
we predict by simulations can indeed be found experi- 
mentally. 
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